Formation spectra of light kaonic nuclei by in-flight (K, N) reactions with a chiral 
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We study theoretically the in-flight (K~ ,N) reactions for the formation of light kaonic nuclear 
systems to get deeper physical insights on the spectra, and to investigate the formation spectra of 
the reaction which will be observed at new facilities like J-PARC. We show the expected spectra 
for the formation of the K~pp, K~ pn, K~nn and K~- 11 B systems which are accessible by the 
(K~,N) experiments. By considering the conversion part of the Green's function, we show the 
missing mass spectra of the (K~,N) reactions in coincident with the particle emissions due to K 
absorption. To calculate the cross sections, we use the so-called Tp approximation to evaluate the 
optical potential. As for the amplitude T, we adopt the chiral unitary amplitude of KN channel 
in vacuum for simplicity. The effects of the p-wave optical potential of E(1385) channel and the 
contributions from K° mixing in 3 He(K~ , n) reaction are also evaluated numerically. We also study 
the behavior of the poles of kaon Green's function in nuclear matter. We conclude that 3 He(if~ , n) 
and 3 He(if~,p) reaction spectra in coincident with the 7rE emission may show the structure in 
the kaon bound region indicating the existence of the unstable kaonic nuclear states. As for the 
12 C(K~ ,p) spectra with the 7rE emission, we may also observe the structure in the bound region, 
however, we need to evaluate the medium effects carefully for larger nuclei. 

PACS numbers: 25.80.Nv, 36.10.Gv, 13.75.Jz, 21.85. +d 
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I. INTRODUCTION 

Kaon-nucleus bound systems, such as kaonic atoms 
and kaonic nuclei, have been studies to obtain informa- 
tion on the kaon properties in the nuclear medium or 
kaon-nucleon interactions at finite density. This infor- 
mation is important as a basic piece of strangeness in 
nuclei. It is also important, for example, for heavy ion 
physics and astro-nuclear physics to determine the con- 
straints on kaon condensation in high-density matter. 

The A'~-nucleus interaction has been studied so far 
based on the experimental data of lightly bound kaonic 
atoms obtained by the X-ray spectroscopy It is 

known that interesting features of kaon-nucleus bound 
systems are tied to the properties of kaons in nuclei which 
are strongly influenced by the change undergone by the 
A(1405) in the nuclear medium, because the A(1405) is 
a resonance state just below the kaon-nucleon threshold. 
In fact, there are studies of kaonic atoms carried out by 
modifying the properties of the A(1405) in the nuclear 
medium @, S 3 ■ These works reproduce the properties 
of specific kaonic atoms reasonably well. The theoret- 
ical studies on kaon-nucleus optical potential based on 
the chiral unitary model were also reported 0, and 
our understanding of the properties of A(1405) has been 
developed @, H, HLE3], recently. It is expected that, in 
near future, one completes the study on the kaon prop- 
erties for lightly bound kaonic atoms accessible by the 
kaonic X-ray technique. The frontier of this subject is 
in higher density region which will be studied by deeply 
bound kaonic atoms [ll[ and kaonic nuclei. 

In recent years, there have been many researches in the 



studies of kaonic nuclear states, which are kaon-nucleus 
bound systems by the strong interaction inside the nu- 
cleus. From the theoretical study of the structure of the 
kaonic nuclear states using ^-nucleus optical potentials 
obtained in the phenomenological approach [l3| and 
the chiral unitary model || , it has turned out that kaonic 
nuclear states have large decay widths of the order of 
several tens of MeV. The kaonic nuclear states for lighter 
HI G3, EH and baryon resonances with two 
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2fj | were also studied theoretically. 



nuclei 
kaons 

Common achievement of the theoretical study on the 
K-pp system E 0, [H EDI is that the K~pp system 
forms a bound state with a rather large width. But the 
binding energy and width do not reach agreement within 
the theoretical calculations yet. In the experimental side, 
an indication of the K~ pp b ound state was observed by 
the FINUDA experiment [21j. However, the critical anal- 
yses of the latest data were reported by Oset, Toki and 
their collaborators [22], [23[ , who claimed that the origins 
of the structure seen in the experimental spectra plj can 
be explained by the well-known processes. The detailed 
analysis was also performed to understand fully the FIN- 
UDA data [Hj]. Thus, the studies of the kaonic nuclear 
structure still have some sources of controversy and we 
need further studies. 

As for the formation and observation of the kaonic nu- 
clear states by direct reactions, in-flight (K, N) reactions 
were proposed and performed by Kishimoto and his col- 
laborators [25|, [2^| . The theoretical calculations of the 
energy spectra of the in-flight (K, N) reaction were per- 
formed in the effective number approach [27l ] and in the 
Green's function method [H, UK, H3- In Ref. (27[, the 
same theoretical technique as Ref. [28[ was used, which 
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was successfully applied for the deeply bound pionic atom 
formation reactions [29|, [3(|. In Ref. [HJ, clear peak 
structure indicating the kaonic bound states has not been 
found in the formation spectra and it has turned out 
that it is difficult to observe clear signals for kaonic nu- 
clear states solely with missing-mass spectroscopy due to 
the large K absorption width. The theoretical predic- 
tion that there are no clear peaks for the K ~- 11 B bound 
system in the formation spectra is consistent with the 
data in Ref. [35| . which observed missing mass spectra 
of the (K~,N) reactions with 1 GeV/c kaon. Based on 
the data of this experiment, Ref. [35[ extracted the po- 
tential depth of the kaon in nuclei by fitting the forma- 
tion spectra in the Green's function method with sim- 
ple background estimation and found that the depth was 
as deep as —160 to —190 MeV. As a new experiment, 
the 3 Hc(A'~,n) reaction for the formation of the K~pp 
system is proposed [361] in the new facility J-PARC. We 
may be able to have new experimental information in 
the coming experiments in near future which will make 
us possible to reach the final conclusion for the proper- 
ties and/or existence of the controversial kaonic nuclear 
states. 

As described above, there are several controversies 
about properties of the light kaonic nuclei. In this ex- 
ploratory level, the purpose of this article is to show the 
formation spectra of various light kaonic nuclei which are 
accessible by the (K, N) reactions, namely K~pp, K~pn, 
K~nn, K°pn, K°nn and K~- ll YS systems, in order to 
compare formation spectra observed in forthcoming ex- 
periments. For this purpose, we study semi-exclusive 
(K~ , N) spectra observed in coincidence detection of the 
final nuclcon with the particles emitted by K absorption 
into nuclei. This is one of the main novelties of this work. 

To calculate the semi-exclusive (K~,N) spectra, we 
exploit the Green's function method with a simple K 
optical potential obtained by a low-density Tp approxi- 
mation. The optical potential is given by the KN scat- 
tering amplitude T and the nuclear density distribution 
p{r). The KN scattering amplitude is evaluated in cou- 
pled channels of KN, 7rE, 7rA, r/Y,, 77A and A'S based 
on the chiral unitary approach and the A(1405) is dy- 
namically generated. This potential automatically in- 
cludes one-body K absorption to the 7rS and irA chan- 
nels. Thus wc can calculate the semi-exclusive processes 
of the kaonic nuclei decays including KN — > 7rE and 7rA. 
We assume a simple two parameter Fermi function for 
the nuclear density p(r), which should be calculated dy- 
namically by few-body treatments though. As we will 
see later, even with such a simple estimation of the op- 
tical potential, the obtained potential produces a similar 
bound state structure with one calculated by an elabo- 
rated few-body formulation. Thus, we expect to obtain 
very similar formation spectra with the few-body formu- 
lation. 

In the present article, we also discuss the effects of 
the p-wave optical potential of the £(1385) channel and 
the contributions from K° mixing in 3 Hc(K~ ,n) reac- 



tion. As for the KN amplitude T, to see the effects of 
the medium modifications, we compare the spectra calcu- 
lated with the in-vacuum amplitude and those obtained 
with the in-medium kaon self-energy evaluated beyond 
the low density approximation based on the chiral uni- 
tary model. 

TABLE I: Possible KNN systems and accessible (K~,N) re- 
actions. Optical potentials used in the calculation of reaction 
are also listed. 

System Formation reaction Optical potential 



K pp 
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K~pn 
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K~nn 


t(K-,p) 


Eq. 




K°pn 


3 Ue(K-,n) 


Eq. 




K°nn 


t(K-,n) 


Eq. 
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In Sec. HIl wc describe the theoretical models for the 
studies of the formation of the KNN systems and the K 
optical potential with the chiral unitary KN amplitude. 
In Sec. IIII) we study the behavior of the poles of kaon 
Green's function in nuclear (proton) matter. Numerical 
results of the formation of the KNN and K~- U B sys- 
tems are presented and discussed in Sec. IIV1 We give 
conclusions of this article in Sec. |Vl 

II. FORMALISM 

To investigate the formation of the light kaonic nu- 
clei such as KNN systems, we apply basically the 
same theoretical framework as our previous publica- 
tions [THHH • We solve the Klein- Gordon equation to ob- 
tain the eigenenergies of the states and apply the Green's 
function method [3l| to calculate the reaction spectra. In 
this section, we explain our theoretical model for study- 
ing light kaonic nuclear systems in detail. 

We consider the (K~,N) reactions as the production 
reaction of light kaonic nuclei with a 600 MeV incident 
energy of K~ in the lab. frame. In this reaction, the in- 
cident kaon kicks out one of the nucleons in the target 
nucleus and the kaon-nucleus system is produced. The 
emitted nucleon is observed in the forward direction to 
make the momentum transfer smaller. Because the inci- 
dent kaon has large momentum, the impulse approxima- 
tion may be good to evaluate the formation spectrum. 
We use the local density approximation to evaluate the 
medium effects on kaon and do not consider the multi 
step process that the kaon hits many nucleons in daugh- 
ter nucleus. 

The energy conservation of this reaction in the lab. 
frame is written as En- + Ma = Ea-i®k- + En, where 
E K - and Ma indicate the energy of the initial kaon and 
the mass of the target nucleus, and Ea-i®k- an( l 
are the energies of the kaon-nucleus systems and the 
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emitted nuclcon. The momentum of the kaon-nuclcus 
system, equivalently the transferred momentum, is in- 
dicated as q. The momentum conservation is written 
as Pk~ = Q + Pn, where Pk- an d Pn indicate the mo- 
menta of the initial kaon and the emitted nucleon, respec- 
tively. In the present reaction, the transferred momen- 
tum is about 200 MeV/c in the lab. frame. We assume 
that the recoil kinetic energy of the kaon-nucleus system, 
T = \q\ 2 /(2Ma-i®k-), is small and negligible. Writing 
the target nucleus mass as Ma = Ma-i + M N — with 
the one nucleon separation energy Sjv, and the kaon en- 
ergy measured from the in-vacuum kaon mass (binding 
energy) as E = Ma-i®k- — Ma-i — M K -, we obtain 
TV = T K - — Sn — E, where T K - is the incident kaon 
kinetic energy, and Tjv the emitted nucleon kinetic en- 
ergy. 

In the impulse approximation, the expected spectra of 

the (K~,N) reactions ( — — ) are evaluated by the 

V dildEN ' 

nuclear response function S(E) and elementary cross sec- 
/ da\ cl ° 
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(1) 



We use the Green's function method [3JJ] to calcu- 
late the response function of the light kaonic nuclei 
in the (K~ , N) reactions. The details of the appli- 
cation of the Green's function method are found in 
Refs. [lj M, EE H S El. The calculation of the 
nuclear response function with the complex potential is 
formulated by Morimatsu and Yazaki [3l| as 



S(E) = --Jm^2 I drdr'T}G{E;r,r')T f , 
71 f 



(2) 



where the summation is taken over all possible final 
states. G(E;r,r') is the Green's function of kaon in- 
teracting in the nucleus and defined as, 



G(E;r,r') = (a\cl> K (r] 



1 



E - H K 



<t>t(r')\a) , (3) 



where a indicates the proton hole state and Hk indi- 
cates the Hamiltonian of the kaon-nucleus system. The 
amplitude tj denotes the transition of the incident par- 
ticle (K) to the nucleon-hole and the outgoing nucleon, 
involving the nucleon-hole wavefunction ipj N and the dis- 
torted waves Xi an d X/i of the projectile and ejectile. By 
taking the appropriate spin sum, the amplitude ry can 
be written as, 

Tf(r) = x}(r)g / 2, m A Y i* i (f)®Tpj N (r)}jMXi{r) , (4) 

with the meson angular wavefunction Y\ R (f) and the spin 
wavefunction £i/2,m a °f tne ejectile. We assume the har- 
monic oscillator wavefunctions for ipj N with the range 
parameter determined by the r.m.s radii of the target 
nuclei. 



The semi-exclusive spectra can be calculated by de- 
composing the response function ^ into the escape and 
conversion parts: S = S csc + S con . This decomposition 
can be done exactly by 

S esc (E) 

= -- I drdr'rj(l + G^vl t )lmG Q (l + V opt G)r f 
n i J 

S con (E) = --Y, I drdr'T}GhmV opt GT f . (5) 
71 f J 

where V op t is the kaon-nucleus optical potential given in 
Hamiltonian. The conversion part is known to express 
the contributions of the kaon absorption to the (K , N) 
spectra [3l|. We can further decompose the conversion 
part of the response function (and clearly the (K, N) 
spectra) into different particle emission channels as, 



S con (E)=Y,SiH(E) , 



(6) 



where j indicates the decay channel. In this work we 
consider j = KN, 7rE, irA, 77A channels. The sepa- 
rated spectrum to each channel can be compared with 
the experimental spectrum observed in coincident mea- 
surements of emitted particles due to the kaon absorp- 
tion. 

To calculate the Green's function G(E;r,r'), we use 
multipolc expansion, 

G(E;r,r')=J2YZ m (n r )Y Lm (n)G L (E;r,r') , (7) 



with the solution G L (E;r,r') of the radial equation, 

. d 2 2d L(L + 1) n _ j 
dr z r dr 



~2vV coul + 2nE)G L (E;r, r') 



(8) 



where \i is the reduced mass of the K and nucleus system 
and V cou \ is the Coulomb potential with a finite nuclear 
size. The strength of the energy dependent potential used 
in this article is also evaluated using this kaon energy. To 
solve the Eq. ([8]), the kaon energy E is regarded as an 
external parameter. The solutions of the radial equation 
can be expressed as (4lj . 

G£(g;ry) = wwjk^) UL{r<)VL{r>) ' 

(9) 

with the regular solution u l (r) at the origin, the outgoing 
solution Vl{t) of the homogeneous equation of Eq. © 
for the fixed values of E and the wronskian A(f,g) = 
fg'—f'g where /' denotes the derivative of / with respect 
to r' . Here, we have introduced the symbol r> (r<) to 
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indicate the larger (smaller) radial coordinate variable of 
r and r' . 

For the values of the nucleon separation energies Sn, 
we use the binding energies of 3 He and t, 7.72 MeV and 
8.48 MeV, for the 3 He(K~,n) and t(K~,p) reactions, 
respectively, since the residual two-nucleon systems are 
not bound, while, for the 3 He(A'~ ,p) reaction, we use 
the proton separation energy of 3 He, which is the energy 
difference of p + d and 3 He systems. For the 12 C(K~ ,p) 
reaction, we use the measured values of the S\/2 and P3/2 
proton separation energies of 12 C given in the Table of 
Isotopes [42j |. 

Since we are interested in the calculation of the for- 
mation spectra of the kaonic nuclei, we use the one-body 
Green function of the in-medium kaon, in which the in- 
teractions of the kaon and nucleus are expressed by the 
optical potential V op t- The optical potential is evalu- 
ated simply by the low density approximation, in which 
the optical potential is given by the forward scattering 
amplitude of the kaon and nucleon in vacuum. For the 
description of the KN scattering amplitude, we exploit 
the chiral unitary approach, in which the KN interac- 
tions are summed up in a non-pcrturbative way in the 
coupled-channels formulation and the A(1405) is dynam- 
ically generated in the KN scattering. Since the iterative 
K interactions with one nucleon in the nucleus is already 
taken into account in the construction of the scattering 
amplitude, we include a correction factor (A — 1)/A for 
the nuclear density to avoid double-counting of the iter- 
ation with the same nucleon. Thus, the optical potential 
is given in the low density approximation by 



(10) 



where p is the nuclear density and T(E) is the KN scat- 
tering amplitude calculated in vacuum. The s-wave KN 
amplitude is taken from Ref. [43|. The importance of 
the ener gy d ependence of the interactions is emphasized 
in Refs. [32|, [44| especially for the imaginary part. We 
include the full e nerg y dependence of the interactions in 
microscopic way |43j both for real and imaginary parts 
in this article. 

The optical potentials between K and NN systems, 
then, can be written as, 



V % NN (r,E) = t-V« NN (E) PNN (r) , (11) 



where pnn indicates the density profile of the two nu- 
cleon systems and V S KNN the s-wave potential strength 



given by the s-wave chiral unitary KN amplitude T s as, 

1 1 

2p2 

1 1 fT s (I = 0)+T s {I = 1) 



V K VP = V K°nn = J_ ^ ( Tg (J = Q) + T# (J = j)) _ ( 12 ) 



V. 



K pn yK°pn 



2p 2 



1 



V s K ~ nn = fT s (I=l) 



(13) 
(14) 



Here, / indicates the isospin of KN, and p the reduced 
mass of the kaon and the NN system. 

The optical potential for the K~- u q system is ob- 
tained in the similar way: 



opt 



2/' 



4 T s (I = 0)+T s {I=l) 5 



11 



11 



Ts(I = 1) 



(15) 



The first and second terms represent the K~p and K~n 
interactions, respectively. Here we have assumed that the 
density distributions for the proton and neutron have the 
same radius and diffuscness parameters. 

In this article, we also study the contribution of the 
p-wave optical potential which includes £(1385) effects. 
We rewrite the Eq. fTl]) with the p-wave potential as, 



KNN 
opt 



{r,E) = Uv s RNN p NN -V-V p RNN PNN V) , (16) 



where V p indicates the p-wave potential strength de- 
fined by the p-wave KN amplitude T p by the same isospin 
combinations as in Eqs. (|12[) - (TT4 |) . T p is defined as the 
first term of Eq. (18) in Ref. [45| including factor 3. 

In our model, we assume the density profile of two 
nucleons as the form of the two parameter Fermi function 
with the radius R and diffuscness a. Here we take R pp = 
Rnn = 1.01 fm and a pp = a nn = 0.50 fm for the pp 



and nn systems, which provide p PP (0) 



fm , and R v 



1.40 fm and a v 



= PnniQ) - 0.117 
0.51 fm for the 



pn system, which is a similar with the deuteron density 
distribution. As we will see later in Sec. IIV1 we will try 
other radius parameters and will find that the calculated 
spectra is insensitive to the change of the parameter. We 
do not specify the isospin of the NN systems here. We 
also use the two parameter Fermi density distribution 
for n B as in Ref. [32| with parameters R = 2.23 fm and 
a = 0.44 fm. 

In the following sections, we compare the formation 
spectra calculated with the optical potential given in 
Eq. pop with those obtained with the in-medium kaon 
self-energy calculated beyond the low density approxi- 
mation in the chiral unitary approach done in Ref. [4(|. 
Since the latter potential has higher-order medium ef- 
fects beyond two-nucleon processes, it is not appropriate 
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for the three-body system. But it may be interesting to 
know the medium effects beyond the low density approx- 
imation and the uncertainties of our theoretical frame- 
work. 

To calculate the conversion part of different particle 
emission channels ((6]), we consider the optical theorem 
satisfied by the chiral unitary amplitude T as, 

3 

where j indicates the intermediate state considered 
to calculate the KN — ► KN amplitude T^ n ^^n 
by the chiral unitary model and means j = 
KN, 7rE, 7rA, 7/A, 77 £, and kS, and o~j is the two-body 
phase space of the intermediate state j. We can interpret 
each term in the right hand side as the contribution of 
the each intermediate channel to the absorptive part of 
K optical potential in the Tp approximation. 

In the present calculations, the final state interaction 
effects to the emitted particles from K absorption are not 
considered, which should be included for more quantita- 
tive results. And the only one-body absorption processes 
are evaluated here and the two-body absorption processes 
such as KNN — > NA are not included in this calculation. 

For the simple estimations of kaon bound states for the 
light kaonic nuclear systems, we solve the Klein-Gordon 
equation with local density approximation assuming two 
nucleon density distribution. And we solved numerically 
as in Re fs. FPU I27I [32| , following the method of Oset and 
Salccdo [43| . Since we adopt the theoretical optical po- 
tential in Eq. with the complex energy dependence, 
we solve the Klein-Gordon equation iteratively to obtain 
selfconsistent solution for complex energy E. 

III. KAON POLE IN NUCLEAR MATTER 

In this section, we study the behavior of the kaon poles 
in the nuclear matter before showing the (K, N) spectra. 
We use the s-wave optical potential in Eq. (fTTj) and ob- 
tain the pole energies by solving the Klein-Gordon equa- 
tion in nuclear matter, 

uj 2 -m 2 K -2m K V s RNN (uj)p = , (18) 

where the ttlk is the kaon mass, the p is the nucleon 
density of the matter and u> is written as u> = E + ttik- 
We neglect the kaon kinetic energy and the p-wave optical 
potential terms by postulating the small momentum of 
the bound kaon. 

The calculated pole positions of K~ in the proton mat- 
ter, where p indicates the proton density and no neutron 
is in the matter, are shown in Fig. [T] in the complex en- 
ergy plane. 

As we can see from the figure, we have only one pole 
at p = at uj = nix corresponding to the free kaon. 
At finite density p > 0, there appear two other poles 
which correspond to the A(1405) resonance of the KN 
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FIG. 1: Trajectories of the poles in the complex uj plane of K~ 
propagator in proton matter obtained by solving the Klein- 
Gordon equation Eq. (JTSJ) for various proton density p. In 
vacuum (p = 0), only single pole appears at the energy of the 
free kaon mass. Crosses correspond to the step size of density 
Sp = 0.25/9o with po = 0.17 fm _,i . The vertical lines indicate 
the threshold energies of neutral KN and 7rE channels. 

system. As pointed out in Ref . Q , both of the two poles 
appeared at p > correspond to the A(1405) resonances. 
As increasing the density p from 0, we find that the kaon 
pole at uj = mx moves to the direction to higher Re uj 
and wider width T = — 21m uj, as shown by the dashed 
line in Fig. [TJ One of the A(1405) pole shown by the 
solid line moves to lower Re uj region with smaller width 
r, which indicates the possible existence of the narrow 
peak structure in the (K, N) spectra in finite nucleus. 
In this calculation, we have only included the one-body 
processes for the kaon absorption, and thus the width 
r become smaller in this energy region because of the 
threshold effects of the 7rE decay. The third pole, which 
starts from uj = (452.2, —66.0) MeV, moves inside the Im 

uj = —65 80 MeV region for p = - 0.5po and never 

comes closer to Re uj axis. 

From the behaviors of the poles in proton matter, we 
may have a schematic picture of the (K, N) spectra. 
Since the kaon pole at w = tuk at p = gets larger 
width r ~ 80 MeV and larger real energy at p = 0.5po, 
this pole is expected to contribute to the quasi-elastic en- 
ergy region in the spectra with large width as a smooth 
slope. The pole shown as the solid line in Fig. [l] which 
is one of the A(1405) poles, has smaller width at larger 
density and, is expected to be seen as peak structures in 
the reaction spectra, which could be considered as can- 
didates of the signals of K~pp states. As seen in Fig. [T] 
the state (pole) is originated from the one of the A(1405) 
poles in the matter and are not connected to free kaon 
at p — » limit. This is one of the interesting findings in 
our framework. In the experimental spectra, this single 
pole may provide several peak structures, since in quan- 
tum mechanics we have in general more than one discrete 
bound level in the finite system for a single pole in the 
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infinite system. The third pole only moves around Im 
w = -65 - -80 MeV (T = 130 - 160 MeV) region, and 
thus, will not make any clear structures in the spectra. 
For further qualitative analyses, one needs to calculate 
the optical potential beyond the low density approxima- 
tion and treat the meson self-energies in a self-consistent 
way as done in Refs. [48| . 

IV. NUMERICAL RESULTS 

A. Kaon optical potential and bound states for 
finite systems 

We show firstly in Fig. [2] the s-wave optical potential 
defined in Eq. (TiTj) to see the strength and the energy 
dependence of the microscopic potential. The contribu- 
tion of the p-wave potential will be discussed later. Since 
we use the microscopic KN scattering amplitude of the 
chiral unitary model, we can consider the potentials of 
the K~-pp, K~-pn, and K~-nn systems simultaneously. 
The energy dependence of each potential appears theo- 
retically in both real and imaginary parts. 

As we can find easily in Fig. [2l the microscopic optical 
potentials have strong energy dependence and channel 
dependence. For the real part of the potential shown in 
the upper panels in Fig. [21 we can see that the K~-pp 
potential is repulsive at the threshold energy correspond- 
ing to the negative sign of the K~p scattering length, 
while it shows the strong attractive nature at lower ener- 
gies due to the existence of the A(1405) resonance. The 
real potential for K~-pn has rather milder energy de- 
pendence than the K~-pp case, and it is almost energy 
independent for the K~-nn system by reflecting there 
are no baryon resonance in the s-wave K~n system in 
this energy region. As for the imaginary potential, only 
one-body absorption processes are included for all cases 
shown in Fig. [21 and thus the strength of the absorption is 
significantly weaker for smaller kaon energies due to the 
smaller phase space of 7r£ decay. As for the p-wave poten- 
tial introduced in Eq. (Tl"6"j) . the strengths of the vv at 
the nuclear center are Re V p K ~ pp (0) ~ 2.5 x 10" 4 [MeV" 1 ] 

and Im V p K ~ pp (0) 0.8 x 10~ 4 [MeV" 1 ] at the kaon 

threshold energy 

We have solved the Klein-Gordon equation with these 
potentials in self-consistent manner for the complex kaon 
energy in finite nuclear systems. The eigenenergies ob- 
tained in the calculation are compiled in Table [ill As the 

TABLE II: Calculated binding energies and widths in unit 
of MeV for the K~ NN systems using the optical potential 
defined in Eq. The widths are shown in parenthesis. 



[MeV] 


K pp 


K pn 


K nn 


Is 


23.0(41.3) 


12.9(36.4) 




2s 


9.0(36.8) 






2p 


11.5(38.5) 


8.6(35.4) 





consequences of the strong complex energy dependence 
of the optical potential, we have found interesting spectra 
of the bound states in K~pp, where the 2s level is sig- 
nificantly close to the 2p level. Because the depth of the 
energy dependent potential becomes so deep at the en- 
ergy of 2p state in complex energy plane and the centrifu- 
gal force is considered in the Klein-Gordon equation, the 
kaon-nucleon subsystem has the component with I =/= 0. 
As we can expect from the potential shape in Fig. [2l 
the number of bound states is different for each system. 
We have also evaluated the effects of the p-wave optical 
potential and found that the variations of the binding en- 
ergies and widths due to the p-wave potential are lesser 
than 10% for the states shown in Table [Til In a U cases 
considered here, the self-consistent solutions in the com- 
plex energy plane have large decay widths, and thus the 
pole position in the complex plane is far from the real 
axis. 

We also notice that since we have assumed the pmn 
profile in this article, we should consider the results listed 
in Table [UJ as qualitative. These binding energies and 
widths could be interesting to interpret the formation 
spectra calculated with the same density distributions 
and interactions. Even in our simple framework, the cal- 
culated results shown in Tabic [TTJ resemble to the results 
of the bound states obtained in three-body calculations 
with a variational approach [l6T |. in which they use a 
KN effective interaction derived by the chiral unitary ap- 
proach [49( | . Since the bound state structures are similar 
with each other, the formation spectra obtained in the 
present calculation could be like the spectra calculated 
with the three-body wavefunction obtained in Ref. [l6T |. 

In Fig. [31 we show the optical potential for the K~- 
pn system evaluated by the K~ selfenergy in the sym- 
metric nuclear matter [46j with medium modifications of 
the chiral unitary amplitude. We see that the real part 
of the potential is attractive even at the threshold by 
the medium effects and has smaller energy dependence 
than that shown in Fig. [2] (b). The imaginary potential, 
which includes the multi-nucleon absorption effects, also 
has milder energy dependence. 



B. Inclusive spectra of KNN system formation 

We show the calculated inclusive 3 He(K~ ,11) spectra 
in Fig. [4] (a) for the formation of the K~pp system with 
the s-wave optical potential defined in Eq. (TiT]) . We find 
that there is a certain bump structure in the subthreshold 
region which could be identified as the signal of the K~pp 
system. Actually we can find the eigen solution of the 
Klein-Gordon equation in this energy region. This pole is 
expected to be connected to A(1405) resonance in p — > 
limit as we have discussed in Sec. lIIII Thus, it is expected 
that this structure can be the signal of the mixture of the 
K~pp and A(1405)p states. It would be interesting to 
study the origins of the structure in the spectra in detail 
as reported recently in Ref. [5fj| for r\ mcsic nucleus case 



7 



40 



> 

£ 

4-1 

Q_ 
O 
> 

2-40 





"a -20 



a. 
o 

> 



40 



-60 



\(a) K-pp 


(b) K-pn 


(c) K"nn 


/ / 
•-' / 
/ 

/ 


— 

■ ■ * ' / 




/ s 
/ / 
/ / 
/ / 
/ / 
^ / 







2 

r [fm] 



2 

r [fm] 



2 

r [fm] 



FIG. 2: The s-wave K"~ optical potentials defined in Eq. fTTJ with the free space chiral KN amplitude T [43| in the Tp 
approximation as a function of the radial coordinate r for the (&)K~ — pp, (b)K~ —pn, and {c)K~ — nn systems, respectively. 
The upper and lower panels show the real and imaginary parts. The solid, dashed, dotted lines indicate the potential strength 
for the real kaon energies Re E = Re u) — ran — MeV, —50 MeV, and —100 MeV with Im E = 0, respectively. 



in connection with the A r *(1535) state in nucleus. 

In Figs. S] (b) and [4] (c), we consider other KNN sys- 
tems accessible by the (K,N) reactions. In Fig. [1] (b), 
we consider the K~pn system formation by 3 He(K~ ,p) 
reactions. Since the optical potential of K~-pn is less at- 
tractive than the K~-pp case in the subthreshold region 
as shown in Fig. [51 the bump structure in Fig. 2] (b) is 
smaller than that in Fig. [4] (a). This subthreshold bump 
for K~ -pn formation, however, still reflects the existence 
of the pole of Klein-Gordon equation in the complex en- 
ergy plane and is interesting to be studied experimen- 
tally. In Fig. |4] (c), we show the expected spectra for the 
formation of the K~nn system by t(K~,p) reaction. As 
we can expect from the weakest K~-nn interaction in all 
systems considered here, we do not find the subthreshold 
structure in the spectra in this case. We only expect to 
observe a possible cusp structure at the threshold which 
could be a reminiscent of bound states. 

The calculated results of the formation spectra for the 
KNN systems in Figs. @] (a)-(c) show the characteristic 
features indicating the differences of the structures and 
optical potentials for the K~pp, K~pn, and K~nn sys- 
tems. It would be interesting that the measurements will 
be performed for all three channels and compared each 
other to see the differences. The signals appearing in the 
inclusive spectra as enhancements in the bound energy 
region, however, are small portions of the whole spectra 
in general and thus, we also consider the semi-exclusive 



spectra later in this section. 



Here we discuss the isospin relation among the KNN 
formation spectra. Since we use the isospin symmetric 
optical potentials as given in Eqs. (fT!?]) and flU]) and we 
do not consider the NN correlations, namely the isospin 
dependence in the NN subsystems, in the present calcu- 
lation, the Green's functions are same for the K~pp and 
K°nn systems, and for the K~pn and K°pn systems, 
respectively. This means that, if there are bound states, 
the spectra of the bound states are equivalent in each 
pair of the systems. But, since the elementary cross sec- 
tion and the distortion factor are different in each system, 
the formation spectra obtained with these potentials are 
quantitatively different in the magnitude. For example, 
we show the calculated spectrum for the K°pn system in 
the 3 He(K~ , ri) reaction in Fig.0](d), which is compared 
with the K~pn formation spectrum in Fig. U (b) calcu- 
lated with the same optical potential strength. These 
spectra have very similar shape qualitatively, but are dif- 
ferent in the absolute values. We can learn the influences 
of the elementary cross section and the distortion factor 
to the formation spectra from these figures. We have also 
confirmed that the spectrum of the K°nn formation in 
the t(K~,n) reaction resembles that of K~pp shown in 
Fig. SI (a). 
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FIG. 3: The K~ — pn s-wave optical potential of the chiral 
unitary approach based on the K~ selfenergy in the symmet- 
ric nuclear matter [46|. as a function of the radial coordinate 
r. The upper and lower panels show the real and imaginary 
parts, respectively. The solid, dashed, and dotted lines in- 
dicate the potential strength for the kaon real energies Re 
E = Re lu - m K = MeV, -50 MeV, and -100 MeV with 
Im E = MeV, respectively. 



Semi-exclusive spectra of KNN system 
formation 



We show in this subsection the semi-exclusive (K, N) 
reactions spectra in coincident with the particle pair 
emissions due to kaon absorption in KNN — > MBN 
one-body processes, where M and B indicate the octet 
mesons and baryons, respectively. We indicate the spec- 
tator nucleon as N in the final state. 

First, we consider the conversion part of the formation 
spectra defined in Eq. (JTJ) with the S con in Eq. ([5]) [3l| . 
which corresponds to the (K, N) spectra accompanied by 
the kaon absorption. The calculated results are shown in 
Fig. [5] for the four kinds of the KNN systems. By com- 
paring these results with the total spectra in Fig. [51 we 
have found that the spectra in the kaon bound region are 
more clearly seen in the conversion spectra. For exam- 
ple, we can compare Fig. [5] (a) with the Fig. [4] (a) for 
the K~pp formation. Since the total spectra in the kaon 
positive energy region include the K~ escaping contri- 
bution, which is removed by considering the conversion, 
the conversion spectrum in Fig.[5]Ja) has prominent bump 
structure in the subthreshold region. Thus, by measuring 



the conversion spectra, we can expect to obtain clearer 
signals in (K, N) spectra. 

Next, we proceed the step further. We calculate semi- 
exclusive spectra by dividing the conversion spectra into 
each particle state as formulated in Eqs. I© and (|17[) . 
The results of the semi-exclusive spectra are shown in 
Fig. [6] for the K~pp system formation. As we can see 
from the figure, by measuring the Especially, the peak 
structures appearing in the subthreshold region in the 
ttT, emission spectra are prominent. This is because the 
K~pp bound states are driven mainly by the strong at- 
traction in the 1 = KN channel with the A(1405) res- 
onance, which decays into the 7r£ state. Thus, our cal- 
culated results indicate that the He(K~,n) reaction in 
coincident with the ttY, emission by K~p — > 7rS one-body 
kaon absorption is one of the best choice to observe the 
subthreshold strength of the spectra to obtain the infor- 
mation on the K~pp state formation. In the K~p and 
K°n emission channels, we have the contributions only 
for the quasi-free energy region, since the KN emissions 
below the threshold are kinematically forbidden. 

In the 3 Ke(K~ , n) reaction, both the K~pp and K°pn 
states are produced and experimentally these states can- 
not be separated out in the observation of the formation 
spectrum. Thus, two contributions should be summed up 
when we compare the calculated spectrum with experi- 
mental observation. In the impulse approximation, the 
spectra for the K~pp and K°pn systems are incoherently 
summed, and the result is shown in Fig. [7] In this figure, 
we plot the spectrum as a function of the emitted neutron 
energy, which is calculated by including the recoil kinetic 
energy of the kaon-nucleus system, since the kaon pro- 
duction thresholds are different in the K~pp and K°pn 
systems. From the figure, we find that the signals of the 
bound energy region in the conversion spectra become 
smaller and less clear by including the effects of K°pn 
formation. However, still there are certain strength in 
this region, which could be observed in the experiments. 

For the comparison of the theoretical calculations of 
the formation spectra with experimental observations of 
the KNN system, it is interesting to consider semi- 
exclusive spectra of two-baryon emissions, which are 
planed to be observed in the experiment performed at J- 
PARC [361 ] . The two-baryon emissions correspond to the 
non-mesonic decays of the KNN system, or equivalently 
the two-nucleon absorptions of K. The non-mesonic de- 
cay of the kaonic nuclei has been recently discussed in 
Rcf. 5l|. The semi-exclusive spectra of the YN emis- 
sions will be discussed elsewhere [531. 



D. Uncertainties of theoretical spectra 

To see the effects of the medium modifications of the 
chiral unitary amplitudes, we show in Fig. [8] the calcu- 
lated 3 Kc(K~ ,p) spectra for the formation of the K~-pn 
systems using the optical potential shown in Fig.[3]bascd 
on the K~ selfenergy evaluated beyond the low density 
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FIG. 4: Calculated spectra for the formation of (a) K~pp in 3 YLe(K~ , n), (b) K~pn in 3 He(K~ ,p), (c) K~nn in t(K~ ,p), and 
(d) K°pn in 3 He(K~ ,n) reactions at T K - — 600 MeV (Pk- = 976 MeV/c) are plotted as a function of the real kaon energy at 
#Jv b = (deg.) for the s-wave chiral unitary optical potential. The horizontal axis E indicates the real kaon energy. Solid line 
shows the total spectra. Dashed, dotted and dot-dashed lines indicate the contributions from kaon s, p, and d partial waves in 
the final state, respectively. The vertical dashed line indicates the kaon production threshold. 



approximation in Ref. (46| for symmetric nuclear matter. 
By comparing Fig. [5] with Fig. 0] (b), we find that the 
subthreshold behavior in Fig. [5] is smoother than Fig. 0] 
(b) and does not show the bump structure because of 
the larger imaginary potential in this energy region. The 
selfenergy used to obtain Fig. [H] includes the higher-order 
medium effects such as multi-nucleon K~ absorption pro- 
cesses in the matter, although the K~-pn system consid- 
ered here, we need to evaluate the medium effects up to 
two-body processes. Thus, we expect that the effects of 
the imaginary potential would not be so strong as seen 
in Fig. 2] and that the medium effects would not change 
the spectra significantly, but the signal could be a little 
unclearer due to the inclusion of two-body absorption. 

We also study the proton density profile dependence 
of the 3 Hc(/v ~, n) total spectra for the K~pp formation, 
fn the previous spectrum calculations shown in Fig. 0] to 
Fig. we have used the density profile given by the Fermi 
distribution with the radius R pp = 1.01 fm and the dif- 
fuseness a pp = 0.50 fm. We also calculate the formation 
spectra of the K~pp systems with the density distribu- 
tions which have the central density with 50% lower and 
50% higher. These density configurations are achieved 
in the Fermi distributions with R pp = 1.58 fm and 0.68 
fm, respectively, which are obtained by adjusting the ra- 
dius parameters with keeping the normalization and the 
diffuscness. The calculated formation spectra with these 
density profiles are shown in Fig. [9] together with the 
spectrum obtained in the previous calculation. As seen 
in the figure, the spectra are insensitive to the 50% varia- 
tion of the central density. We have carried out the same 



study for the K~nn formation spectra and confirmed 
that the spectra are also insensitive to the change of the 
density distribution. Having also performed the study of 
the density profile dependence of the formation spectra 
with a different K optical potential, which is given in 
Ref. [53l | , we have obtained the similar formation spectra 
within the variation of the density distributions. 

In Fig. 1101 we show the contribution from p-wave opti- 
cal potential, which takes account of the £(1385) res- 
onance effects in KN channel. We have used the p- 
wave scattering amplitude calculated in Ref. (45| in chiral 
unitary model and used the optical potential defined in 
Eq. (|16|) . We found that the contributions from the p- 
wave optical potential is tiny for the total spectra in our 
present theoretical framework as in Fig. [TD] This results 
is consistent with that in Ref. f54|, where the p-wave 
potential effects are very small for kaonic atoms. The ef- 
fects of the p-wave potential and / or £(1385) resonances 
are also reported in Refs. [l5|, [l6|, [55[ in the studies of the 
structure of the KNN systems, where these effects are 
not negligible for quantitative discussions. 

Finally, we would like to add a few comments on the 
uncertainties of our results. In our framework, only one 
body operator is considered as shown in Eq. ((4]), and the 
momentum transfer in the (K, N) reaction is shared by 
the high momentum components of the wavefunctions of 
a nucleon in the target in the initial state and a kaon 
in the final state as in Ref. (34[. The many body re- 
action processes, which we have not evaluated here, ex- 
ist in reality and have certain contributions to the re- 
action spectra. We think that the contributions from 
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FIG. 5: Conversion parts of the calculated formation spectra of (a)K~pp in 3 He(K~ , n), (b)K~pn in 3 He(7f~,p), (c)K~nn in 
t(K~ ,p), and (d)K°pn in 3 He(K~ ,n) reactions at = 600 MeV (Pg = 976 MeV/c) are plotted as a function of the real kaon 
energy at #5v b — (deg.) for the s-wave chiral unitary optical potential. The horizontal axis E indicates the real kaon energy. 
Solid lines show the total conversion spectra. Dashed, dotted and dot-dashed lines indicate the contributions from kaon s, p, 
and d partial waves in the final state, respectively. The vertical dashed line indicates the kaon production threshold. 



the many body processes are structureless and mainly 
appear in the lower emitted proton energy region in gen- 
eral, which corresponds to the quasi-free kaon produc- 
tion region. Hence, we think the calculated results in the 
quasi-free kaon production region can be affected by the 
many body processes j35j . Several processes may have 
finite contributions to the spectrum in the bound kaon 
energy region as backgrounds as discussed in Ref. [44| . 
which will be evaluated quantitatively in future works. 

E. (K~ ,p) spectrum for 12 C target 

Finally, we show the calculated spectra of the 
12 C(K-,p) n B®K- reaction at T K - = 600 MeV with 
the chiral unitary amplitude. In Fig. [TTJ we show the to- 
tal and the conversion spectra calculated with the chiral 
amplitude in free space, together with the total spectra 
with the amplitude at finite density for comparison [32| . 
As we can see from the Fig. [Til the spectra around the 
threshold becomes smooth and structure-less due to the 
medium effects to the chiral unitary amplitude. This is 
the same tendency which we have observed for the K~pn 
systems in Fig. [4] (b) and Fig. [8] The medium effects 
which include the multi-nucleon processes for kaon ab- 
sorption make it more difficult to observe clear structures 
in the total spectrum. 

We show in Fig. [12] the semi-exclusive spectra for the 
12 C(K~ ,p) n ~B®K~ reaction calculated with the chiral 
amplitude in free space. These figures show the (K~,p) 
spectra with the meson and baryon emissions due to kaon 



absorption as explained for the K~pp system case in 
Fig. [H In the 12 C target case, we have two different fea- 
tures from the K~pp system, which are the contributions 
from p 3 /2 proton hole state in 11 B and the kaon absorp- 
tion by neutron. In Fig. I12[ we show the contributions 
from different KN isospin and proton-hole states sepa- 
rately. The contributions of kaon absorption by neutron 
are shown in three lower panels. We find that the quali- 
tative features of the semi-exclusive spectra are resemble 
to the K~pp case and we find again that the 7r£ emission 
channel from the initial K~p subsystem have large com- 
ponents in the bound energy region. We also find that 
the 7r~A, 7r~£° and 7r°£~ emission channels from K~n 
initial subsystem have certain strength in subthreshold 
kaon energies. Though we do not include the final state 
interaction for the emitted meson and baryon after kaon 
absorption, it would be interesting to explore the semi- 
exclusive spectra even for larger nuclei such as 12 C since 
we have possibilities to observe clearer structure in the 
spectra and to obtain information on subthreshold kaon 
properties at finite density, which could not be observed 
in total spectra alone. 



V. CONCLUSION 

We have made systematic studies for the formation 
spectra of the KNN and K'-^B systems which are ac- 
cessible by the (K, N) reactions. We have adopted the 
theoretical KN amplitude obtained by the chiral unitary 
model and formulated the optical potential within the so- 
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FIG. 6: Semi-exclusive spectra for the formation of the K pp system in 3 He(_ftT , n) reactions at T K - = 600 MeV (Pk- = 976 
MeV/c) at 6£ ab = (deg.) for the s-wave chiral unitary optical potential. As indicated in the figure, each figure corresponds 
to the different meson(M)-baryon(B) emission channel after one-body kaon absorption process K~ N — » MB. Solid line shows 
the total spectrum for each semi-exclusive channel. Dashed and dotted lines indicate the contributions of KN isospin= and 
1 optical potential, respectively. The vertical dashed line indicates the kaon production threshold. 
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FIG. 7: Calculated results of 3 Re(K~,n) reaction spectra at T K - = 600 MeV (P K - = 976 MeV/c) including both K~pp and 
K°pn in the final states are shown in (a) total spectra and (b) conversion part at #^ ab = (deg.) for the s-wave chiral unitary 
optical potential. Dashed and dotted lines indicate the contributions from K~pp formation and K°pn formation, respectively. 
Solid lines are the sum of the both contributions. 



called Tp approximation to calculate the formation spec- 
tra systematically. This theoretical optical potential has 
the complex energy dependence in contrast to the phe- 
nomenological potentials for which we need to assume (or 
neglect) the energy dependence. Although we have used 
the Green's function with the density dependent opti- 
cal potential for the KNN systems rather than treating 
them as few-body systems, this simplification makes the 
spectrum calculation much simpler. 

In order to know the basic features of the light kaonic 
nuclear systems, we solve the Klein-Gordon equation 



with the theoretical potential sclfconsistently for kaon in 
proton matter and for KNN bound states with assuming 
the two nucleon density distributions. From the behav- 
iors of poles in the proton matter, we found that the pole 
in the kaon bound energy region is connected to one of the 
poles of A(1405) in low density limit and not connected to 
free kaon in vacuum. The solutions of the Klein-Gordon 
equation of the KNN systems show clearly the differ- 
ences of the structure of the bound spectra for K~pp, 
K~pn, and K~nn states due to the different characters 
of the interactions. 
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FIG. 8: Same as Fig. [4] (b), except for the K~ — pn optical 
potential used in the calculation. The K~ selfe nerg y II in the 
symmetric nuclear matter with medium effects [4g| is used as 
the K~-pn optical potential. 
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FIG. 9: Calculated results of 3 He(K , n) spectra at T K - — 
600 MeV (P K - = 976 MeV/c) for the formation of K~pp 
state at #J; ab = (deg.) with the s-wave chiral unitary optical 
potential. Solid line shows the same result as that in Fig. [4] 
(a). Dashed and dotted lines indicate the results with the pp 
distributions with R pp — 1.58 fm and 0.68 fm, which have 
50% lower and 50% higher central densities, respectively. See 
details in text. 



As for the formation spectra of the light kaonic nuclei, 
we have made detail analyses along to the line described 
in Sec. HI We would like to emphasis that the semi- 
exclusive spectra for (K, N) reaction coincident with the 
particle emissions due to kaon absorption are important 
to obtain the experimental information on the subthresh- 
old kaon properties at finite density. Especially, in 7rS 
emission channels we expect to obtain subthreshold kaon 



properties and clearer indications of kaon bound states 
as enhancements in the spectra. Namely it is possible to 
obtain the information on the existence of bound states 
by knowing the strength of the subthreshold spectrum. 
On the other hand, it is difficult to know the binding 
energies and widths quantitatively from the spectra in 
general for the systems with large width. 

Finally, we mention the effects of the mixture of 
the K°pn formation process to K~pp formation in the 
3 Hc(K~ ,n) reaction. This effect should be evaluated 
since it is included experimental spectra as a kind of back- 
ground. We have reported the calculated results in our 
theoretical framework in this article. 

We believe that our theoretical results are interesting 
and stimulating for both theorists and experimentalists 
in this field and help to develop the research of the kaon 
physics. 
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FIG. 10: Same as Fig. |4](a), except for the K~ — pp optical 
potential used in the calculation. Dashed line indicates the 
result calculated with the s-wave chiral unitary optical po- 
tential (Eq. (HU), and solid line indicates that with s- and 
p-wave optical potential ( Eq. fl!6j l). The p-wave chiral am- 
plitude calculated in Ref. [45| is used to evaluate the optical 
potential. 
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FIG. 12: As indicated in the figure, each figure corresponds to the different meson (M)-baryon(i?) emission channel after 
one-body kaon absorption process K~ N — > MB. Due to the isospin symmetry, the spectrum with 7r°£~ emission channel is 
the same as that of 7r~E° channel. Upper and middle six panels show the emission channels from K~ +p and lower three panels 
the emission channels from K~ + n, respectively. Thick solid line show the total spectrum for each semi-exclusive channel. 
Thick and thin dashed lines indicate the contributions of KN isospin=0 with sj/2 and p 3 / 2 proton hole states, and thick and 
thin dotted lines indicate those of KN isospin=l with s 1 / 2 and p 3 / 2 proton hole states, respectively. There are no isospin=0 
contribution in the K~ + n amplitude (lower three panels). 
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